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Abstract 

We establish strict upper limits for the Casimir interaction between 
multilayered structures of arbitrary dielectric or diamagnetic materials. 
We discuss the appearance of different power laws due to frequency- 
dependent material constants. Simple analytical expressions are in good 
agreement with numerical calculations based on Lifshitz theory. We dis- 
cuss the improvements required for current (meta) materials to achieve a 
repulsive Casimir force. 

Dated: 20 Oct 2005, Europhysics Letters, in press. 

PACS. 42.50.Pq - Cavity quantum electrodynamics; 42.50.Lc - Quan- 
tum fluctuations, quantum noise; 78.67.-n - Optical properties of low- 
dimensional, mesoscopic, and nanoscale materials and structures 

1 Introduction 

The optical properties of materials that show both a dielectric and magnetic re- 
sponse, have recently attracted much attention (see ^ for a review). A number 
of striking phenomena hke perfect lensing and a reversed Doppler effect have 
been predicted, and experimenters have begun to explore the large parameter 
space of structural units that can be assembled into artificial materials. Break- 
throughs have been reported on the way towards designed susceptibitilies in 
the near-infrared and visible spectral range j2| Ej ■ Quantum electrodynamics 
in meta materials has recently been explored with particular emphasis on left- 
handed or negative- index materials 01 |S] . We discuss here to what extent the 
Casimir interaction between two meta material plates can be manipulated by 
engineering their magneto-dielectric response. Strict limits for the Casimir in- 
teraction are proven that apply to all causal and linear materials, including both 
bulk and multilayer structures. We illustrate these results by computations of 
the Casimir pressure, considering materials with frequency-dispersive response 
functions like those encountered in effective medium theories. We derive power 
law exponents and prefactors and find that a strongly modified Casimir interac- 
tion is possible in a range of distances around the resonance wavelength of the 
response functions. We give estimates for the required temperature range and 
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structure size: it is not unreasonable to expect that improvements in fabrication 
and detection will allow for experimental observations. 

One of the most striking changes to the Casimir interaction is a cross over 
to repulsion. This has been predicted previously for idealized magnetodielectric 
materials [SlIZIEl and objects suspended in a liquid [HlEIl. In the latter case, 
repulsion has been observed experimentally with colloidal particles ^Jj and is 
also used in a recent proposal for measuring Casimir torques ■ Casimir repul- 
sion between mirrors separated by vacuum requires a strong magnetic response 
[3 El that hardly occurs in conventional ferromagnets |14l I15| . Indeed, to ma- 
nipulate the Casimir force in the micrometer range and below, where it can be 
conveniently measured, the key challenge is to achieve a magnetic susceptibility 
at high frequencies, approaching the visible range. Now, there is a well-known 
argument due to Landau, Lifshitz, and Pitaevskii that ij-{lo) = 1 in the visible 
This objection, however, only applies to materials whose magnetization 
is of atomic origin, where the magnetic susceptibility is Xm ~ {v/cY' ^ 1. An 
array of split ring resonators with sub-wavelength size typically gives, on the 
contrary, Xm ~ {'-^MY'I ^ li where i? is the resonance frequency and / the 
filling factor |2I E] ■ As we illustrate below, artificial materials that are struc- 
tured on the sub-micron scale are promising candidates for a strongly modified 
Casimir interaction. 



2 Lifshitz theory 

For two perfectly conducting plates held at zero temperature and separated by 
a distance d, Casimir derived a force per unit area given by Fc = 7r^?ic/(240 d"*) 
[T!?j . We use the convention that Fc > corresponds to attraction. For linear 
media with complex, frequency-dependent material parameters, the force can 
be computed from Lifshitz theory |2U| . This expression has been re-derived, 
for plates of arbitrary material and for multilayer mirrors, using different meth- 
ods EHESESlEllEHlEZIEHlEniEniEIlES- At finite temperature, it can be 
written in the form 

where the sum is over the imaginary Matsubara frequencies ujn — i^n = 2TTinkBT/h 
(the n = term being multiplied by 1/2), and k is related to the wave vector 
component perpendicular to the mirrors, kz = (w^/c^ — — kyY^^ = in. The 
i^\a (A = TE, TM, a = 1, 2) are the reflection coefficients at mirror a for elec- 
tromagnetic waves with polarization A |22l 1281 . For homogeneous, thick plates, 
they are given by 



'"TM = = (2) 

e(iCn)c^ + v/a(e(iCn)Mi60 - 1) + ^^2^2 

(exchange e and n for tte)- The zeros of Dx = e'^'^'^/{rxir\2) — 1 at real 
frequencies define the eigenmodes of the cavity formed by the two mirrors. 
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3 Strict limits 



To derive upper and lower limits for , we use that the Kramers-Kronig rela- 
tions |16| imply real and positive material functions at imaginary frequencies, 
^(i^) ^ 1: provided the material is passive (non-negative absorption Ime(u;) > 
0). As a consequence, the Fresnel formulas (0) imply — 1 < rxa < 1, and we 
find 

111 

~ e^'^'i + 1 ~ Ih ~ e2'''i - 1 ^ ' 

with the stronger inequalities < 1/Dx < l/{e'^^'^ — 1) holding for identical 
plates. In the latter case, the Casimir force is hence necessarily attractive. The 
inequalities (O saturate for a perfectly conducting mirror facing a perfectly 
permeable one (si = oo, /i2 = oo, say), and for identical, perfectly reflecting 
mirrors, respectively. The resulting forces at zero temperature are |2()l IHj 

T = 0: -Ifc <Fl< Fc. (4) 

o 

In the high-temperature limit, we get similarly 24 — |i^T < < Ft = 
C{3)kBT / (SttcP) by keeping in Eq.lQJ only the n = term in the sum. 

Consider now a mirror made from layers of arbitrary passive materials. Re- 
flection coefficients for such a system can be obtained recursively. For a layer 
'6' separating a medium 'a' from a substrate 'c', for example, 

fabc — -, 9it-t,,, Wi 

where Vab (^6c) describes the reflection from the interface ab (be), respectively, 
and w is the layer thickness |33[ |2| • If the substrate c is a multilayer sys- 
tem itself, Tfcc is the corresponding reflection coefficient. For the imaginary 
frequencies occurring in the Lifshitz expression Q, the wavevector in the layer 
is purely imaginary, kh ~ iki,, and single- interface coefficients are real [Eq.Q]. 
From Eq.jSJ, they remain real for multilayer mirrors. In addition, the mapping 
fab fabc is & conformal one, and if r^c e^^fib-i" jg j-gg^j g^j^(\ g [—1,1], the interval 
[— 1, 1] is mapped onto itself. For multilayer mirrors, we thus obtain again the 
inequalities — 1 < < 1. This generalizes the limits of Refs. "311122 that are 
obtained only for layered dielectric mirrors, using transfer matrices. 



4 Casimir interaction between metamaterials 

To illustrate these generally valid results, we focus on mcta materials described 
by effective medium theory El E] • We adopt Lorcntz-Drude formulas for e 
and /i 

MiO = i + ^, ^^(^0 = 1 + ^. (6) 

Regarding the permeability, we have taken the limit of weak absorption and 
computed /i(iO in terms of Im/i(a;) using the Kramers-Kronig relations. This 
is necessary to ensure high-frequency transparency of the medium. We denote in 
the following by a typical resonance or plasma frequency occurring in Eqs.®. 
The corresponding wavelength, A = 27rc/r2, provides a convenient distance scale. 
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Note that a (magnetic) resonance wavelength as short as ~ 3 /im has already 
been achieved with material nanofabrication The key advantage of meta 
materials is that their electric and magnetic 'plasma frequencies' Via and Qa 
are fairly large as well: a value of Qa ~ "da^/J < (c/a)-\/7 is typical for a 
split-ring resonator array with period a and filling factor / jl7| . This property 
is also necessary, of course, to achieve a left-handed medium (eQ,(w), ^a{^) < 
for some real frequencies). The magnetic plasma frequencies occurring in 
conventional ferromagnets are much smaller |14| , and the impact on the Casimir 
interaction is weak, as reported recently 

In the plots shown below, the Casimir pressure is normalized to hfl/cP [see 
after Eq. 121)]. In order of magnitude, this corresponds to 10^pNmm^^/(A//xm)'* 
at a distance d = A/2. This can be measured with sensitive torsion balances [201 

or cantilevers |38[I39| . We plot in FigQlthe result of a numerical integration 
of Eq.(^, the curves corresponding to different material pairings. One sees 
that in all cases, the force satisfies the limits (@J that exclude the shaded areas. 
We observe that materials with negative index of refraction around show a 
strongly reduced attraction (Fig^Jb)). This can be attributed to the reduced 
mirror reflectivity due to impedance matching. Casimir repulsion is achieved 
for some distances between mirrors made from different materials (Figmc,d)). 
At short distance, i.e. d <^ A/27r, even these pairings show attraction with a 
power law 1/d'^. Coating one mirror with a magnetic layer (Fig^c)), there 
is a sign change around the layer thickness w: for A/2tt ^ d <^ w, the layer 
behaves like a thick plate, and its material parameters lead to repulsion. The 
layer can be ignored for w d, and one recovers the attraction between the 
(identical) substrates. This is consistent with asymptotic analysis based on the 
reflection coefficient as we outline below. Detailed calculations show that 
a large resonance frequency is not sufficient to achieve repulsion, the oscillator 
strength of the resonances (proportional to ili and Q2) must be large enough so 
that hfli,hQ2 ^ Taax{kBT,hc/d). As the temperature is raised, the distance 
range where repulsion is observed disappears, see Fig|5] One then finds a 1/d^ 
power law at large distance as well. 

The different regimes of Fig^can be understood from an asymptotic analysis 
of Eq.QJ. At short distance {d ^ A/27r), the integral is dominated by a region 
in the K-^-plane where the Fresnel coefficients Q take the nonretarded forms 
''TM ^(e) = (e— f)/(e-t-l)>0 assuming that e > I and similarly tte 
R{fJ,) > unless /i = 1. Proceeding like Lifshitz [2010111; yields to leading order 
a power law — cs/d^ with a positive Hamaker constant given by 

^3 = ^ ^'{Li3[i?(£l(60)^(£2(en))] +Li3[i?(A*l(60)^(A*2(en))]}, (7) 

n=0 

where Li„(z) = J^'kLi ^'^ / ■ must be noted that for the special case of 
homogeneous plates, this asymptotic expression actually provides another, much 
stricter, upper limit to the Casimir force, since rxMa < Ri^a), fTEa < RifJ-a) 
and Li3(z) is a monotonous function (see Fig^. In order of magnitude, C3 ^ HQ, 
at low temperatures {ksT <^ hn). Compared to ideal mirrors, dispersive plates 
thus show a much weaker Casimir interaction that is in general attractive (Fig^ 
and Ref. 35 ). At larger distances, A/2tt ^ d <C At = hc/ksT, the Casimir 
force follows a 1/d^ power law, and repulsion is found provided one of the 
materials is dominantly magnetic. Here, the non-dispersive results of Ref. [7] are 
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Figure 1: Casimir force between planar mirrors of different dispersive materials 
at zero temperature. Left: plot of the force per unit area, normalized to hQ,/cP 
vs. the plate separation scaled to A = 2-kc/VI. Here, J7 is a typical plasma or 
resonance frequency. Shaded areas: excluded by the inequalities (0}. Right: 
material response functions along the real and imaginary frequency axis. 
Dashed lines: short-distance asymptotics c^/d^ with coefficient {T)), providing 
another upper limit for homogeneous mirrors. Curve (a) (inset): two identical 
non-magnetic Drude metals (in Eq.®, f^a = Vt, uja — 0, Ha{0 = 1, ck = 1,2). 
Curve (b): two identical left-handed meta materials with overlapping dielectric 
and magnetic resonances {ila = 0.3 il, uja ~ ^, Oa — 0.3 51, — fl). Curve 
(c): two identical non-magnetic Drude metals one of which is coated with a 
left-handed meta material (Drude metals with fli = 3fi, = 0, /ii(C) = 1; 
left-handed coating with thickness w = 10 x 27rc/51 and material parameters 
(dominantly magnetic response) = 0.1 £7, UJ2 = 51, O2 = 0.3 51, '02 = 51)- 
Curve (d): two meta materials, one purely dielectric, the other mainly magnetic 
(51i ^3n,LUi^n, ^i(^) = 1, fl2 = 0.151, LU2 = 57, 62 = 0.3 51, t?2 = 57). 

recovered. Finally, for d 3> Ay, the leading order force is the term n = in the 
sum , again an attractive 1/d'^ law, with a coefficient given by an expression 
similar to (O, but involving the static material constants, see |7]. 

The impact of temperature is illustrated in Fig|21 at high temperature, 
fc^T ^ hn, the second Matsubara frequency ^1 falls already into the mirrors' 
transparency zone, and the 1/d^ power law is valid at all distances. As T — > 0, 
the intermediate repulsive zone appears in the range A/27r <C d ^ At/2t:. 
A good agreement with the analytical 1/d^ asymptotics is found outside this 
zone, as shown by the dashed lines. For the resonance wavelength A = 3 fim 
mentioned above, cooling to a temperature T « 0.1 hn/kB ~ 50K is required 
to 'open up' the repulsive window. This temperature increases, of course, with 
materials whose response extends to higher frequencies. 

Finally, we would like to illustrate the kind of peculiar asymptotics that 
becomes possible with carefully matched material parameters. This follows 
Ref.|25 that computes the Van der Waals force on a water film coated on 
both sides by lipid membranes, finding a weak dependence on the ultraviolet 
frequency range because both materials have a similar electron density. Consider 
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Figure 2: Casimir force between two different meta materials as the temperature 
is raised. The material parameters are identical to Line (d) of FigQ] The tem- 
perature takes the values fc^T = (a)0.3, (6)0.1, (c)0.03, {d)Q Tifl. Same scaling as 
in Fig^ Dashed lines correspond to the short- and long-distance asymptotics 
discussed in the text. 

thus a liquid-filled gap with a similar electron density as medium 2 so that 
£o = £27 and a permeability /zq = ^1 = 1 matched to medium 1. For simplicity, 
we assume that these equalities hold at all frequencies. In this case, we can 
show that the force is repulsive at all distances, even at finite temperature. 
Indeed, both contributions in Eq.Q vanish, and the leading order term for high 
temperatures also vanishes. The high-temperature limit is given by the n = 1 
term in This gives a distance dependence proportional to exp(— 47rd/AT) 
similar to what has been observed in some experiments with colloids (mentioned 
in 01]). As the temperature is lowered, this exponential regime still applies for 
d ^ At. If ksT ^ hU, the regime of Ref. f?" exists at intermediate 

distances A/27r <C d <C At/2tt. The short distance regime sets in for d <C A, 
and an analysis similar to the one leading to Eq.Q gives (T = 0) 

00 

n ^ if;(2„.f-/|r(3-2„.2„f./.) (-£i-|<«^) (8) 

n-l g 

with T{k,z) = f^dtt'^^^e"'^ . At short distance, the sum is dominated by the 
first term, so that to leading order, we get a repulsive power law Fl — —ci/d 
(Figin^a)). In order of magnitude, ci ~ hft^ /c^ and therefore again —Fl ^ Fc, 
with a cross over occurring around A/27r (see Fig|3Ja)). Due to our assumption 
of a perfect matching £2 = eoi this kind of behaviour seems quite remote from 
experimental reality. As shown in Fig|3Jb-d), a slight mismatch between the 
dielectric functions of liquid and plate leads back to an attractive force, first at 
short distances, then suppressing the repulsive window altogether. 

5 Conclusion 

We have generalized strict upper and lower limits for the Casimir force. We 
have shown that a strongly modified Casimir force can occur between dispersive 
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Figure 3: Casimir force between two different meta materials as the residual 
dielectric response of mirror 2 is changed. Same scaling as in Fig^ T = 0. 
Dashed lines: short-distance asymptotics The dielectric response for mir- 
ror 2 is of Lorentz-Drude form with fixed cutoff frequency and variable oscillator 
strength. Material functions as in Fig^d), except £2(0) = 1 (a), 1.01 (b), 1.03 
(c), 1.1 (d). 



and absorbing mirrors with a sufRciently large magnetic susceptibility, extend- 
ing results restricted to non-dispersive materials The most promising way 
to achieve this repulsion seems the use of meta materials engineered at scales 
between the nanometer and the micron because they provide a fairly large mag- 
netic oscillator strength. Our results arc intrinsically limited to distances d 3> a 
by our use of effective medium theory. Sufficiently small structures and suffi- 
ciently low temperatures then ensure that in the range a <C rf <C AT/27r, the 
Casimir interaction can be strongly altered: even if repulsion cannot be achieved 
in a first step, we expect a significant reduction of the Casimir attraction at dis- 
tances of a few microns (Fig^. 
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